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A detailed description of the baryon direct Urea processes A: n —> p + e + t'e, B: A ^ p + e + 

C: —>■ A + e + J/e related to the neutron star cooling is given in the relativistic mean field 

approximation. The contributions of the reactions B and C on the neutrino luminosity are calculated 
by means of the relativistic expressions of the neutrino energy losses. Our results show that the total 
neutrino luminosities of the reactions A, B, C within the mass range 1.603-2.067M0 (1.515-1. 840Mq 
for TMl model) for GMl model are larger than the corresponding values for neutron stars in npep 
matter. Although the hyperon direct Urea processes B and C reduce the neutrino emissivity of the 
reaction A, it illustrates the reactions B and C still make the total neutrino luminosity enhancement 
in the above mentioned areas. Furthermore, when we only consider the proton superfluidity 
in neutron star cooling, we find that although the neutrino emissivity of the reactions A and B is 
suppressed with the appearance of ^5o proton superfluidity, the total contribution of the reactions 
A, B, C can still quicken a massive neutron star cooling. These results could be used to help prove 
appearing hyperons in PSR J1614-2230 and J0348+0432 from neutron star cooling perspective. 

PACS numbers: 21.60.-n, 26.60.-c, 26.60.Dd, 24.10.Jv, 13.75.Cs 


I. INTRODUCTION 

Neutron star (NS) constitutes one of the best astro- 
physical laboratories for studying dense matter physics. 
It arises at the end of life of a 8-20 Mq massive stars and 
forms in the aftermath of the core collapse supernovae ex¬ 
plosion. A newly born NS is very hot with temperature 
as high as lO^^-lO^^K, but rapidly cools to a temperature 
of less than 10^°K within minutes. The cooling process 
of a NS is dominated by a combination of surface pho¬ 
ton emission and interior neutrino emission. The latter is 
responsible for about 10^-10® years until the interior tem¬ 
perature reaches 10®K. It is generally known that photon 
luminosity is obviously lower than neutrino luminosity, 
meaning that the thermal radiation from a NS surface 
reflects the intensity of interior neutrino emission [I|, 
While neutrino emisision depends strongly on the compo¬ 
sition of super dense matter in NSs. It is well known that 
NSs cores are dense enough to allow for emerging exotic 
matter with the strangeness quantum number through 
weak equilibrium, such as A, S°, E+, E“, 5°, 5“ hy¬ 
perons, referred as nphe/i matter, except for the con¬ 
ventional nucleons and leptons (npe/r matter) 1H3. It 
means that all the possible baryon neutrino emission pro¬ 
cesses would happen during the neutrino cooling stage 
Among them, the most powerful enhancement 
of neutrino emission is provided by the nucleon direct 
Urea processes (NDUP), secondarily is the hyperon direct 
Urea processes (HDUP) (2ll - [^ . Besides, the degrees of 
freedom of hyperons tend to soften the equation of state 
(EOS) calculated in the relativistic mean held (RMF) 
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model based on SU(6) spin-flavor symmetry (quark model 
for the vector meson-hyperon coupling constants), then 
reduce the maximum mass of NS to about 1.6-1.7Mq 
[29l - [3^ . However, Demorest et al. in 2010[3^ indicated 
that the binary millisecond pulsar PSR J1614-2230 ex¬ 
panded the maximum observational mass from 1.67 ± 
O.O2M0 to 1.97 ± O.O4M0 using the Shapiro delay mea¬ 
surements from radio timing observations. Antoniadis et 
al. in 2013(33 observed another massive neutron star 
PSR J0348-I-0432, whose mass is 2.01 ± O.OAMq. It is 
clear that the inclusion of hyperons in such heavy NS 
cores are difficult to explain by SU(6) spin-flavor symme¬ 
try in RMF model. And for this reason, the SU(3) fla¬ 
vor symmetry is widely applied to RMF model. Because 
it changes the strength of the isoscalar, vector-meson(a; 
and (j)) couplings to the octet states, which can sustain 
a NS with mass of (1.8-2.1)M0 even if hyperons exist 
in NS core [ssl - l^ . Furthermore, baryons in NS interior 
can become the superfluid state related to the generation 
of BB Cooper pairs under attractive interaction. The 
baryon superfluidity(SF) could suppress considerably the 
NDUP, HDUP and thus affect the cooling rate of NS re¬ 
markably [3, [U, ( 43 . As we all know, the neutrons in the 
crust and protons, hyperons in the core undergo Cooper 
pair in state, while neutrons in the core can pair in 
^P 2 state. 


This paper is arranged as follows. In section 2, we 
make a brief review for RMF, NS cooling theories and 
the gap equation for the ^Sq proton SF. The numerical 
results are discussed in section 3. Finally, we summarize 
our conclusions in section 4. 
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TABLE I. The parameter sets GMl and TMl. The relations, Qct^n = gpA = 0i assumed. We take niu, = 783MeV, 
nip = 770MeV, mjv = 938MeV. For the GMl and TMl models, mo- = 550MeV and 511.198MeV, respectively. 


gujh guiT. guiB g<t,A g^s g^s 

GMl SU(6) 7.073 7.073 3.537 5.002 5.002 10.003 

GMl SU(3) 8.149 8.149 6.038 -6.253 -6.253 -9.004 

TMl SU(6) 8.409 8.409 4.205 5.945 5.945 11.891 

TMl SU(3) 9.689 9.689 7.180 -7.435 -7.435 -10.706 

TABLE 11. The other coupling constants for hyperons. The relations, gpN = = gpH, are assumed. 


II. THE DENSITY EQUATIONS 
A. RMF Theory 


above Lagrangian into the Euler-Lagrange equations 


dL ^ dL 


( 2 ) 


In this calculation, we adopt RMF model to describe 
NS matter. The constituents of NSs fall into two cat¬ 
egories: npe/r and nphe/r matter. The strong interac¬ 
tion between baryons is mediated by the exchange of 
isoscalar scalar and vector mesons cr, w, isovector vec¬ 
tor meson p. The two additional strange mesons are also 
included, namely isoscalar scalar cr* and vector </> mesons 
[13, [13, • The total Lagrangian is given by 

L = - QaBCr - ga-Ba>) ( 1 ) 

B 

-gpBip.T ■ p^ - gcpBjfj.uj'" - g<j>BjtJ.(l)'^]tpB 
+ ^{9pCrd^a - m^cr^) -h ^{dya*d''a* - 

- \p^'’Ppy + 

+ \rnlp^p^ + -f 

Here — dp^ujy Rp,y — dpPy dyp^ and Pp,y — 

dp,4>y — dy(l)pt denote the field tensors of w, p and tp mesons, 
respectively. The sum on B and 1 runs over the octet 
baryons and leptons, namely, n, p. A, E+, S°, 

S“, e, p. ipB, 'Pi and ms, mi are the baryon, lepton Dirac 
helds and masses, respectively. 7 ^ is the Dirac matrice. 
The meson helds are replaced by their expectation values 
at the mean held level. Now we are going to plug the 


The equations of motion for each baryon and meson helds 
can be obtained in RMF approximation 


- ga,Bjoi^° - PpBIoTopI 
-50B7O<)>°)'*/'B = 0, 

(3) 

gaBPSB = mla + acr^ -|- ba^, 

B 

(4) 

guiBPB = mlujo + caWp, 

B 

(5) 

''^PpBPbIsb = m^pPos, 

B 

(6) 

'^^ga-BPSB = 

B 

(7) 

'^gpBPB = m^^ipo- 

(8) 


B 


Here Jb and I^b express the baryon spin and isospin 
projections, respectively. m*g is the baryon effective mass 

rriB = mB - gaBcr - ga-Ba-, (9) 


The scalar density psB and baryon density pB are given 

by 
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c 

h 
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cos 9c 

1 
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-7^(EaD/3) 

sin 9c 

^3/2 

7372(E - D/3) 

cos 9c 

1 

F-D 

sin 9c 

-1 

D-F 

cos 9 c 

0 


cos 9 c 

72 

V2F 

sin 9c 

7172 

{F + D)/V2 

sin 9c 

1 

Fad 


TABLE III. The constants of the baryon direct Urea processes. We take sin 0^=0.231A0.003, F=0.477A0.012, D=0.756A0.011. 


The hadron phase should meet the local charge neutrality 
and beta-equilibrium conditions. The former is given by 

Pp + Pi:+ = Pe- + Ph- + Pe + Pp- (12) 

The latter is imposed by the baryon chemical potential, 
which is a linear combination of /i„ and /ie, 

PB = Pn-qBPe,Pe= Pp (13) 

where qb is the baryon electric charge(in unit of e). 

We can solve the Eqs. (3)-(13) self-consistently at a 
given baryon density pB- 

The total energy density and pressure of NS matter 
are 



1 ,_, 1 rPFi , - 

+ + pi dpi 

i -^0 



Eqs. (14) and (15) as inputs, we can obtain the mass- 
radius relation by sol ving the Tolman-Oppenheimer- 
Volkoff(TOV) equation |44]. 


dP{r) 

dr 


[P{r) A £{r)][M{r) A 47rr^P(r)] 


r{r — 2M{r)) 
dM (r) 


dr 


= 47rr^e(r). 


(16) 


We adopt two successful RMF parameter sets to de¬ 
scribe NS matter, GMl and TMl, as listed in Table I 
[ 3 ^. These parameters have been determined by fitting 
to some ground state properties of nuclear matter. As 
for the couplings of the isoscalar vector mesons uj and 
(j) to baryons, we adopt two relations: SU(6) spin-flavor 
symmetry based on the naive quark model and general 
SU(3) flavor symmetry, as listed in Table II (40l| . 


B. NS cooling theory 


The baryon direct Urea processes consist of two suc¬ 
cessive reactions, beta decay and capture, are listed in 
Table III [H. 


31 — B 2 A e A B 2 A e —>■ Bi A (1'^) 

Here Bi and B2 represent baryons. Due to the EOSs of 
NSs matter are derived by RMF model, so the neutrino 
energy losses must be consistent with the used relativistic 
EOSs. In the free relativistic gas, the energy and momen¬ 
tum conservations require a large effective mass differece 
of Bi and B2, ra*g^ lOOMeV, which is unlikely to 

appear in the reactions A, D and G. The reason is that 
the effective masses of hyperons with the same species 
but the different isospins are same in Eq.(9). Therefore, 
in the relativistic regime, the energy conservation should 
be assured by considering the potential energy difference 
of Bi and 32- The neutrino emissivity can be given by 
the Fermi Golden Rule 

Qo = 2 I [H - /2)(1 - /3)|My,|fl8) 

X { 2 t:)'^S{Ei - E2-£3- £ 4 )^(Pi - P2 - Pa - P4), 

where Pj, £j express the momentum and kinetic energy 
of particle species j (j = 1,2,3 and 4 refer to Bi, B 2 , e 
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and Ve), respectively, fj is the Fermi-Dirac distribution 
functions of baryons and electrons, 

1 


Jb = 


(19) 


exp((£;_B - fJ-B)/T) + 1’ 

f =_i_. 

exp((£:3 - He)/T) + 1 

The delta functions 6 {Ei — E 2 — — £ 4 ) and 

'^(Pi — P 2 ~ Ps — P 4 ) describe the energy and momen¬ 
tum conservation. E 12 = £ 1,2 + C^i .2 is the energy of 
baryons. Ui ^2 is the potential energy of baryons, which 
can be obtained in Section A and has the following form 

Ub = 9uiB^0 + 9pBhBPo + g<t,B(t>o- ( 20 ) 

Namely, 

Un Up = 9pNP0^ (21) 

U-B- —Uy,o = —gpY,Po, 

Ub.- - Ubo = -gpBPo- 

is the squared matrix element of the baryon 
direct Urea processes summed over spins of initial and 
final particles 

= 32GlC^[(gf - ff)M*M;{P4P3) + {91 (22) 

-fl?{P4P2){P3Pl) + (51 + flf{P4Pl)iP3P2)l 

where Pj = {ej,Pj). Gp = 1.436 x 10“"^® erg cm^ is 
the weak-coupling constant, /i, gi and C are the vector, 
axial-vector constants and Cabibbo angle which are given 
in Table 111. 

The relativistic expression of the energy loss Q per unit 
volume and time in NS matter is expressed as [461 147| 

4^77r 

Q = (23) 

X + £F 2 )Pf 3 ~ ~ ^F 2 ){Pfi — PF 2 )) 

+ 2glg,emlm*2 + {fl + gl){g.e{2£F^eF2 - w^m^) 

+ £FiP3 — 2 ^Pfi ~PF 2 PPF3)i^Fi +£j’2))}) 

In this expression, pfi, PF 2 and pf 3 are the Fermi mo¬ 
menta of baryons and leptons. Ef^ and £f 2 are the ki¬ 
netic energy of baryon at the Fermi surface. 0 = 1 if the 
Fermi momenta pfi, PF 2 , PF 3 satisfy the triangle con¬ 
dition and 0 = 0 otherwise. The situation of muons is 
similar to that of electrons. 

The cooling equation based on the approximation of 
isothermal interior is. 


— —{Li, + Lr)- 


(24) 


Here Ln and L^- are the total neutrino and photon lumi¬ 
nosities, respectively. Gy is the total thermal capacity of 
NS matter. They are 


L^= / =47ri?^cr(10r)3e 




C, 


i;0 


= f {Ge + Gp + Y,CBo)dV, 

'j D 


(25) 


where a is the Stefan-Boltzmann constant, e* = 
a/1 — 2 mjr is the gravitational redshift. e^‘ is the value 
of at the stellar surface (r=R). The semiempirical ex¬ 
pression Ts = (ior)®/3 expresses the relation between 
interior temperature T and surface temperature Tg. 


C. SF of protons 


The key quantity in determining the onset of ^S'q pro¬ 
ton SF is the gap function A(p), 


A(p) 


1 

dTT^ 


p-pp np,P)Mp )_ 

yiVT+W) 


(26) 


where £{p) = E{p) — E{pFp)- E{p) is the single-particle 
energy of protons with momentum p 


E{p) = ^Jp^ + m*^ + gyjpuj + g^pfj) + gpphpP- (27) 

V{p,p ) is the pp potential matrix element. In this work, 
we use the Reid soft core(RSC) potential for the pp 
potential , as an example to demonstrate the influ¬ 

ence of hyperons on the ^So proton pairing energy gaps. 
The critical temperature Typ of proton SF is given by 
the pairing gap A(fe) at zero temperature approximation. 


Typ = 0.57 A{k). (28) 

As a result, the neutrino emissivity and thermal ca¬ 
pacity can be written as 


Q = QoRb, Gb = GbqRcbo- (29) 


Here Rb and Rcbq are the SF reduction factors of the 
neutrino emissivity and thermal capacity, respectively. 
For the ^S'o proton SF, the reduction factors Rp and Rcpo 
are 


0.0163exp(^i^^^^) ^ 3.149exp(^i^^^^) 

~ ^ T ’ ^Cpo - / T N2..q 

' Tep ' ' Tep ' 

(30) 

According to the discussion of the RMF approach above, 
we can obtain the EOS, mass-radius relations and neu¬ 
trino emissivities of the reactions A, B, C as well as Fermi 
momentum and single particle energy of protons, then 
the pairing gap and critical temperature of proton 
SF and speed of the NS cooling can be obtained. 


III. RESULTS AND DISCUSSION 

In this section, we give three cases in Eq.(I) for RME 
theory: (i) the non-strange cr, w, p mesons are included 
in SU(6) spin-flavor symmetry; (ii) the cr, w, p mesons in¬ 
cluding strange mesons cr* and cj) are considered in SU(6) 
spin-flavor symmetry; (iii) cr, w, p, a* and tp mesons are 
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FIG. 1. EOSs including hyperons in NS matter. 


taken into account in SU(3) flavor symmetry. We mainly 
study the effects of the degrees of freedom of hyperons 
and reactions B, C on the EOS, neutrino emissivity, neu¬ 
trino luminosity, energy gap of proton SF and NS 
cooling. Then we compare our results with PSR J1614- 
2230 and J0348-I-0432, whose measured masses are used 
as reference values. 

Fig. 1 shows the EOSs in three cases. Fig. 2 shows the 
mass-radius relations of NSs by solving the TOV equa¬ 
tion. The softest and hardest EOSs are obtained by cases 
(i) and (iii), respectively. Though the coupling g^]\f for 
case (iii) is smaller than the corresponding value for case 

(i) as shown in Table 1, the total repulsive force is at¬ 
tributed not only to cj meson but also to (j) meson. As 
seen in Fig. 1 and 2, though we consider the contribu¬ 
tion of the strange mesons cr* and (j) on the EOS in case 

(ii) , the coupling gtj,N=0. It means that (j) meson only 
couples to hyperons and makes the EOS be not enough 
stiff. So the hardest EOS is obtained only through the cj) 
meson in case (iii). From case (i) to (iii), the maximum 
mass of NS (the corresponding center density) sequently 
increases from 1.820 (0.771), 1.863 (0.817) to 2.141M0 
(0.871) for the GMl model, 1.686 (0.673), 1.729 (0.754) 
to 2.O38M0 (0.848) for the TMl model, respectively (Fig. 





FIG. 2. Mass-radius relations including hyperons. 


2). Namely, the EOS in SU(3) flavor symmetry could be 
consistent with the observed values of PSR J1614-2230 
and J0348-I-0432 when hyperons appear in NS core. Fig. 
3 depicts the neutrino emissivities of the reactions A, B, 
C in nphe/r matter for the three cases. In order to make 
the effects of hyperons and reactions B, C on A more in¬ 
tuitive, Fig. 4 depicts the total neutrino emissivity of the 
reactions A, B, C in npe^ and nphe/r matter for case (iii). 
As can be seen from Figs. 3 and 4, the neutrino emis¬ 
sivity of the reaction A has a tendency to decrease with 
increasing of the baryon density pb due to the presence of 
hyperons in NS matter decreases the nucleon and lepton 
fractions according to the requirement of the charge neu¬ 
trality and P equilibium conditions(Eqs.(12) and (13)). 
Also, the neutrino emissivities of the reactions B and C 
are obviously less than that of A, it is due to that they 
have smaller matrix elements in Eq.(22). The strongest 
neutrino emissivity of the reaction A or B is in case (iii), 
while the weakest one is observed in case (i). For the 
reaction C, the neutrino emissivity in case (iii) is less 
than the corresponding values in cases (i) and (ii) firstly 
and then increases slowly, equals or exceeds the values in 
cases (i) and (ii). While the reactions D-F would never 
happen within stable NSs, because their threshold densi- 
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FIG. 3. Neutrino emissivities of the reactions A, B and C as 
a function of the baryon density pB in nphe/r matter. 


ties are larger than the center densities of maximum mass 
NSs. As shown in Figs. 2 and 3, the mass ranges of the 
reactions B and C in case (iii) are 1.671 — 2.141Mq and 
1.888 - 2.141M0 for the GMl model, 1.579 - 2.O38M0 
and 1.849 — 2.O38M0 for TMl model, respectively. 

Next we mainly discuss the effects of the degrees of 
freedom of hyperons and reactions B, C on the total neu¬ 
trino emissivity, neutrino luminosity, energy gap of ^S'o 
proton SF and NS cooling in case (iii). Fig. 5 gives 
the radial distributions of the total neutrino emissivities 
of the reactions A, B, C for the GMl model with 1.98, 
2.00, 2.10, 2 .I 2 M 0 stars in case (iii). The radial distribu¬ 
tions of the total neutrino emissivity for the same mass 
stars in npe/r and nphe/x matter are almost unchanged 
when r is relatively large (Part I). However, the reactions 
B and C happen in succession as the radius decreases 
(Part II and III), leading to the radial distributions of 
the total neutrino emissivities in nphe/r matter are sig¬ 
nificantly larger than the values in npe/i matter. While 
the radius regions of the increasing total neutrino emis¬ 
sivities shrink continually with the increasing NS mass. 
The situation of the TMl model is like the above in GMl 
model. Fig. 6 shows the total neutrino luminosity as a 



FIG. 4. Total neutrino emissivities of the reactions A-F as 
a function of the baryon density ps. The solid line is the 
neutrino emissivity of the reaction A in npep matter. The 
dashed line is the total neutrino emissivity of the reactions A- 
F in nphep matter. The dotted line is the neutrino emissivity 
of the reaction A in nphep matter. 


function of NS mass m in case(iii). As seen from Fig. 
6, whether hyperons are included or not, the neutrino 
luminosity increases firstly and then decreases with the 
increasing of the baryon density ps- Once the mass of 
NS reaches a value, one luminosity corresponds to two 
different NSs. And the total neutrino luminosities of re¬ 
actions A, B, G within the mass range 1.603 — 2.O67M0 
and 1.515 — 1.84OM0 will be larger than the values in 
npe^ matter for the GMl and TMl models, respectively. 
Fig.7 shows the ^S'o proton SF critical temperatures as 
a function of the baryon density ps for case(iii) in npe/i 
matter (solid lines) and nphe/i matter (dashed lines), re¬ 
spectively. One can find that whether or not the NS core 
appears hyperons, the proton SF critical tempera¬ 
ture first increases, reaches a peak and then decreases 
to zero as the baryon density pB increases. While when 
the NS core appears hyperons, the ^S'o proton SF critical 
temperatures are first below and then above the corre¬ 
sponding values in npe/i matter within the density ranges 
of pb = 0.0 — 0.454 fm“^(pB = 0.0 — 0.418 fm“^ for the 
TMl model) and ps > 0.454 fm“^(pB > 0.418 fm“^ for 
the TMl model) for the GMl model, respectively. It is 
because the total contribution of the Fermi momentum 
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r(km) 


FIG. 5. Radial distributions of the total neutrino emissivities 
with different mass NSs in npe/r matter(solid lines) and nphe^ 
matter (dotted lines) in the GMl model 


and single-particle energy of protons (see Eqs.(ll) and 
(27)) results in Tqp change. Besides, the density range 
of the proton SF is widened due to the inclusion of 
hyperons. The range of proton SF can achieve cover¬ 
age or partial coverage in the cores of NSs, which is highly 
relevant to the reactions A, B. That is, the presence of 
hyperons affects not only the reaction A and proton 
SF critical temperature, but also the NS cooling. In Figs. 
8 and 9, the cooling curves of NSs are calculated through 
solving the cooling Fq.(24) by assuming isothermal NS 
cores for case(iii). Observational data of 8 isolated NSs 
whose effective surface temperatures have been measured 
or constrained are listed as compared with the theoret¬ 
ical cooling curves As you can see from Figs. 

8 and 9, the cooling curves of moderate mass NSs can 
explain the observational data, while the cooling curves 
for massive NSs due to the low surface temperature are 
difficult to explain the existing observational data. And 
the cooling curves decrease smoothly with the increasing 
NS mass, meaning that the cooling of massive NSs under¬ 
goes faster neutrino cooling whether the reactions B and 
C appear in NSs. In addition, the same mass stars with 



m 


FIG. 6. Total neutrino luminosity of the reactions A, B, C as 
a function of the NS mass m. The solid line is the neutrino 
luminosity of the reaction A in npe/r matter. The dashed line 
is the total neutrino luminosity of the reactions A, B, C in 
nphep matter. The dotted line is the neutrino luminosity of 
the reaction A in nphe/r matter. 


the reactions B, C(dashed lines) are colder than that in 
npe/i matter(solid lines). Taking case (iii)in GMl model 
as an example, we can see that the neutrino luminosity of 
1.70, 1.95, 2.03 Mq NSs in nphe^ matter are greater than 
the corresponding values in npe/i matter which is because 
the three NSs within the mass range 1.671 — 2.O67M0, 
see Fig. 6 for details. While the proton SF criti¬ 
cal temperatures Tcp of 1.70, 1.95, 2.03 M 0 in nphe/i 
matter is lower than the value in npe/i matter(see Figs. 
2 and 7 for details), in the case we studied here. As a 
result, the reactions A, B are suppressed in advance in 
nphe/i matter. However, the reaction C is not affected 
by the proton SF. Therefore, although the neutrino 
emissivities of the reactions A and B are suppressed with 
the presence of the proton SF, the total contribution 
of reactions A, B, C can still speed up a massive NS cool¬ 
ing. For the TMl model, the situation is similar to the 
described in GMl model. The cooling curves presented 
here are estimates in order to make the effects of the ^S'o 
proton SF on the reactions A and B be more clearly, in 
particular, the effects of NS crust and hyperons SF are 
not considered. 
































FIG. 7. The critical temperature Tcp of the ^So proton SF as 
a function of the NS mass m in wpefj, matter(solid lines) and 
nphe/r matter(dashed lines) 


IV. CONCLUSION 

We have studied the effects of the degrees of freedom of 
hyperons and reactions B, C on the reaction A in NS mat¬ 
ter using the two popular RMF parameter sets, GMl and 
TMl, respectively. Firstly, we used the SU(3) flavor sym¬ 
metry to obtain the relatively stiff EOS for supporting 
the observed massive PSR J1614-2230 and J0348-I-0432. 
Secondly, the total neutrino luminosity of the reactions 
A, B, C are calculated in npe/i and nphe/i matter, re¬ 
spectively. We found that the presence of the reactions 
B and C caused of the total neutrino luminosity higher 
than the corresponding values in npe/r matter within the 
mass range 1.603 — 2.O67M0 for the GMl model and 
1.515 — 1.84OM0 for the TMl model. Finally, our main 


purpose is to test the effects of the ^S'o proton SF on 
the reactions A, B by comparing cooling curves with ob¬ 
served data. We will analyze the effects of other baryon 
SF on the corresponding baryon direct Urea processes in 
future work. Our results showed that the cooling rate of 
the same mass of two NSs with the reactions B and G are 
obviously faster than that without the reactions B and 



FIG. 8. Observational data(error bars) on surface tempera¬ 
tures of 8 NSs as compared with theoretical cooling curves 
obtained by the GMl model for proton SF from Fig.7. The 
solid lines correspond to npe/r matter, the dashed lines corre¬ 
spond to nphe^ matter with masses(from top to bottom)!.5, 
1.7, 1.95 and 2.03Mq, respectively. 

G. These features maybe can help to prove the presence 
of hyperons in massive NSs cores. 
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